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Abstract. Recent precise observations of solar global parameters are used to 
calibrate an upgraded solar model which takes into account magnetic fields 
in the solar interior. Historical data about sunspot numbers (from 1500 to 
the present) and solar radius changes (between 1715 and 1979) are used to 
compute solar variability on years to centuries timescales. The results show 
that although the 11 year variability of the total irradiance is of the order of 
0.1%, additional, longer lived changes of the order of 0.1% may have occurred 
in the past centuries. These could, for example, account for the occurrence of 
climate excursions such as little ice ages. 



Introduction 

Because of the steady increase of the concentration of 
greenhouse gases in the Earth's atmosphere due to an- 
thropogenic activity, global warming is potentially the 
most serious environmental threat of our time. In order 
to assess the magnitude of the threat, we must under- 
stand the sensitivity of the climate system to a variety of 
external driving processes. In particular, the magnitude 
of the solar variability component of climate change, po- 
tentially very significant, is still relatively uncertain. 

One way of assessing the effects of a given driver of 
climate is by trying to fit the observations (e.g. repro- 
ducing the air temperature variations of the last 150 
years [Jones et at, 1986; Hansen and Lebedeff, 1988] or 
sea surface temperature [Reid, 1991]) with, and with- 
out the driver in question operating. For the case of 
the Sun, we can assume that it does not change, and 
adjust all the other parameters to reproduce the tem- 
perature record. Subsequently, a given model of the 
solar variability is assumed, and then all the param- 
eters are re-adjusted to reproduce the same temper- 
ature record. Obviously, if the solar variability con- 
tributes to the observed trend, the model sensivisity 
to the other processes that were assumed to produce 
the trend on their own must decrease. Similarly, if the 
solar component varies in antiphase with the observed 
trend, the sensitivity to the other processes must be in- 
creased. Andronova and Schlesinger carried out such 
a calculation [Andronova and Schlesinger, 2000]. They 
used two models for the solar variability. One, con- 
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structed by Lean et al. [1995], henceforth denoted by 
L, was based on sunspot areas and (when available) lo- 
cations, He 1083 nm emission, group sunspot numbers, 
and Ca emission from the Sun and solar analogues. The 
other model, developed by Hoyt and Schatten [1993], 
(henceforth called HS), was based on the fraction of 
sunspot area occupied by the penumbra, equatorial ro- 
tation rate, length and decay rate of the solar cycle, 
and the mean level of solar activity. Upon including so- 
lar variations, the model sensitivity to greenhouse gases 
changed drastically, decreasing by 53 percent for the HS 
solar variations, and by 44 percent for the L solar vari- 
ation. Obviously, if either model provides an approxi- 
mately accurate representation of solar variability over 
the last 150 years, the solar input to climate variability 
is very important. 

Our Solar Variability Model 

The problem with the models of solar variability re- 
ferred to above is that they are based purely on empiri- 
cal correlations which are not supported independently 
either by theory, or by direct, quantitative observations. 
As a consequence, we do not know how credible these 
extrapolations are. In fact, although the explanation 
that the short term (days to months) variability of the 
total irradiance is largely due to active regions (spots 
and faculae) is beyond question, the precise origin of 
the 11 year modulation is still uncertain. The most 
popular explanation is that the modulation is similarly 
caused by contributions from active regions plus the 
magnetic network, which are known to vary with the 
activity cycle. In this explanation, the photosphcric 
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Figure 1. Comparison between the measured (solid 
curve) and calculated (dashed curve) solar irradiance 
variations. 



background, and everything below the photosphere, re- 
main unchanged. Alternatively, Li and Sofia Li and 
Sofia [2000] propose a different explanation in which the 
1 1 year modulation of the irradiance is due to structural 
adjustments of the solar interior produced by variations 
of an internal magnetic field located at about 0.96i?Q. 
These adjustments modify the energy flow from the in- 
terior to the surface, and affect the entire luminosity of 
the Sun. This would imply that both the structure of 
the Sun and the photospheric temperature should vary 
over the activity cycle, and various observations of the 
solar photospheric temperature [Gray and Livingston, 
1997a; , b], solar oscillation frequencies [Libbrecht and 
Woodard, 1990] and solar radius [Emilio et at, 2000; 
Antia et al, 2000] support this explanation. 

A variable magnetic field in the solar interior adds 
a variable magnetic component to pressure and inter- 
nal energy not considered in standard models. In addi- 
tion, the field modifies the transfer of convective energy. 
As a consequence, the structure of the solar interior 
changes (albeit slightly) from the field-free state, and 
so do all global parameters, such as radius and lumi- 
nosity (hence, effective temperature). The formulation 
of this model upgrade is developed in Lydon and Sofia 
Lydon and Sofia [1995], and updated in Li and Sofia Li 
and Sofia [2000] . Model calculations show that all con- 
ventional (nonmagnetic) global parameters are sensitive 
to the location, intensity, distribution and orientation 
of the solar interior magnetic field, and this enables us 
to use the observed cyclic variations of solar radius, lu- 
minosity and effective temperature to determine cyclic 
variations of those magnetic properties. 
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Figure 2. Comparison between the measured (solid 
curve) and calculated (dashed curves) solar photo- 
spheric temperature variations. 

It is well known that the even splitting coefficients 
vary with solar activity cycle [Libbrecht and Woodard, 
1990; Bhatnagar et al, 1999; Howe ey al, 1999]. Since 
magnetic fields contribute to the even splitting coeffi- 
cients, and the field causes the activity as it buoys to 
the surface, we can relate the maximum magnetic field 
B m to the sunspot number Rz in the following way: 

B m = B {A+[l + log(l + R z )} 5 }. (1) 

The magnetic field inferred from the solar oscillation 
frequencies observed around 1996 was equal to about 
20 kG, and it peaked at r « 0.96i?Q [Antia et al, 
2000]. To match the even splitting coefficients, the 
magnetic field distribution in the solar interior is as- 
sumed to have a gaussian-like profile, and to be purely 
toroidal. Since the location of the field is fixed at around 
r = 0.96i?Q by more accurate helioseimology, the more 
uncertain observation of the solar radius variation is not 
needed. Consequently, B , A, and the profile width can 
be determined by fit for the simultaneously-observed 
cyclic variations of solar irradiance [Wilson and Hud- 
son, 1991; Frohlich and Lean, 1998] and effective tem- 
perature [Gray and Livingston, 1997a; , b]. 

Since 1978, several space-borne radiometers have 
monitored solar irradiance and provided us with daily 
irradiance data of unprecedented precision [Frohlich 
and Lean, 1998]. In our model calculations we use one 
year as the timestep, so we use yearly averages of the 
total solar irradiance. The standard variances range 
from 0.1 to 0.5 Wm" 2 . The measured irradiance varies 
almost in phase with the solar activity cycle with a rel- 
ative amplitude of about 0.1% (solid curve in Fig. 1). 
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Gray and Livingston [1997a, b] measured solar effective 
temperature variations over the 1978-1992 interval by 
monitoring the neutral carbon A5380 line in the solar 
flux spectrum. They found that the temperature varies 
about 1.5 ± 0.2 K during a cycle (solid curve in Fig. 2, 
also nearly in phase with the cycle. Since this line has 
a high excitation potential (7.68 eV), it forms deep in 
the photosphere. As a result, the effective temperature 
measured in this way is actually a photospheric tem- 
perature, thus it is free of the effects of changing ac- 
tive regions and network. The best fit for the observed 
cyclic variations of solar irradiance and temperature is 
obtained by setting B = 90 G, A = 190, and O.O6i? 
for the magnetic field profile width. Our model fits the 
irradiance observation well (xl — 4), and the temper- 
ature a little less well {\\ = 43). To obtain the best 
fit we minimize \ 2 = Xi + X%- This best fit model is 
also consistent with helioseismology [Antia et al, 2000]. 
For example, using Eq. (1) with B = 90 G, A = 190 
and R z = 8.6 in 1996, we find B m w 20 kG. Of course, 
it peaks at r = 0.96i?Q as required by helioseismol- 
ogy [Antia et al, 2000]. The resulting W parameter is 
almost constant and equal to about 2 x 10~ 5 . 

Extrapolation of the Activity Cycle 
Timescale Variability 

After calibrating our model by using the recent pre- 
cise data for solar irradiance during the period from 
1978 to 1999, and for effective temperature during the 
period from 1978 to 1992, we can extrapolate the so- 
lar irradiance back during the period when the annual 
sunspot numbers are available. 

Surely the best-known features of the Sun are the 
sunspots and the regular cycle of the solar activity, 
which waxes and wanes with a period of about 11 years. 
Naked-eye observations of sunspots have been made in 
China since before A.D. 300, and telescopic observa- 
tions in Europe have been made since A.D. 1610. The 
discovery of the sunspot cycle by Schwabe in 1843 led 
Wolf in 1847 to 1851 to the concept of universal sunspot 
numbers, which in 1868 he extended backward annu- 
ally to 1700. At Zurich, Wolf calculated magnitudes of 
each cycle from 1749, and the central dates or epochs of 
maxima and minima back to about 1610. Early aurorae 
or Northern lights provide clues to the dating of early 
sunspot cycles. The printed catalogs and additional 
"Spectrum of Time" material [Schove, 1948, 1967; Eddy, 
1976; Schove, 1961] have made it possible to extend 
Wolf's tables [Wolf, 1851] back into earlier centuries. 

On scrutinizing the historical sunspot data, Schove 
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Figure 3. Solar variability in the past five centuries 
(dotted curve corresponds to the Zurich sunspot num- 
ber Rz) and in the past four centuries (solid curve conr- 
responds to the group sunspot number Rg)- kG stands 
for kilo-Gauss, L for total solar luminosity, T for solar 
effective temperature, R for solar radius. 



Schove [1983] obtained an annual sunspot number table 
starting from A.D. 1500, as shown in the top panel of 
Fig 3 (dashed line). The gap from 1645 to 1715 corre- 
sponds to the so-called Maunder minimum. The second 
panel from the top in Fig. 3 shows the maximum mag- 
netic field strength calculated by Eq. (1) with the best 
fit parameters. Obviously, the magnetic field strength 
equals to about 17 kG even if the sunspot number is zero 
(e.g. in 1810), since helioseismology requires 20 kG for 
1996. The last three panels show the cyclic variations 
of solar luminosity, radius and effective temperature, re- 
spectively. We also use the more precise group sunspot 
numbers [Hoyt and Schatten, 1998] to recalculate solar 
variability in the past, as shown by the solid curves in 
Fig. 3. Obviously, both results have the same trend. 
Because luminosity changes result in solar irradiance 
changes and the relationship is linear, the (relative) so- 
lar irradiance change is the same as the (relative) lumi- 
nosity change. Detailed tabulations of these data will 
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be provided on request. 

Additional Variability Mechanisms 

As we can see from Fig. 3, the maximum variabil- 
ity of the solar radius is about 2 x 1CP 5 , or 0.02 arc s. 
Although this variation is in agreement with the most 
recent determination of the cycle radius variations ob- 
tained from the MDI experiment on SOHO [Emilio et 
al, 2000], it is much smaller than the radius changes de- 
termined from historical data over the last 2 centuries. 
In our view, the most reliable historical data sets from 
which solar radius changes can be determined are the 
duration of total eclipses measured near the edges of 
totality. From them, changes of the order of 0.5 arc s 
have been detected. In particular, a change of 0.34 arc s 
between 1715 and 1979 (Dunham et al 1980), a change 
of 0.5 arc s between 1925 and 1979 [Dunham et al, 
1980], and no change between 1979 and 1976 [Sofia et 
al, 1983], were detected. If such changes are real, what 
could cause them? What are the corresponding solar 
irradiance changes? 

If we use the magnetic field location required to pro- 
duce the 11 year cycle variability, we find that it is im- 
possible to produce a 0.5 arc s radius variation even 
if we apply an unreasonablely strong magnetic field. 
However, our model shows that the deeper the loca- 
tion of a magnetic field and the more intense the mag- 
netic field, the larger the resulting radius change. The 
physical cause is obvious: the stronger the magnetic 
field, the larger the contribution of magnetic pressure 
to the fixed total pressure required by hydrostatic equi- 
librium. Consequently, the gas pressure decreases ac- 
cordingly and the gas will expand to increase the solar 
radius. However, gas expansion is limited by the gas 
density. Therefore, we compute the magnetic field re- 
quired to produce the detected radius change between 
1715 and 1979 as a function of mass depth, as displayed 
in the top panel in Fig. 4. It is well known that a 
strong magnetic field will cause a change of location of 
the boundary between the convective and the radiative 
region [Lydon and Sofia, 1995]. The second panel from 
the top in this figure shows how the convection bound- 
ary Rcz varies with the applied magnetic field (solid 
curve), and how the location of the maximal magnetic 
field, Rb, varies with the mass depth (dashed curve). 
The shadowed region indicates the half-width of the re- 
quired magnetic field. Of particular relevance are the 
values corresponding to the base of the convection zone, 
as indicated by the dot-dashed line in this figure, since 
all conventional dynamo models locate the process pre- 



cisely at that depth. There, the magnetic field required 
to cause a 0."34 change of the solar radius is 1.3 million 
G, and the resulting luminosity variation is 0.12 percent 
(the third panel of Fig. 4), which is almost equally due 
to the variation of effective temperature (the bottom 
panel) and radius, since the radius variation contributes 
2 x Ami? = 0.07%. These values are interesting for 
producing significant climate change if the solar vari- 
ations are sufficiently long lasting, and for not grossly 
contradicting what we know about the Sun, excepting 
a value for the magnetic field that is larger than we are 
comfortable with. 

Results and Discussion 

According to the picture presented in this paper, the 
solar variations of years to centuries are due to varia- 
tions of magnetic fields at two different depths in the 
solar interior. The 11 year cycle would arise from mag- 
netic fields of the order of 20-40 kG, located at or above 
O.96i?0. At the same time, there is a stronger field lo- 
cated at just below the base of the convection zone that 
undergoes changes on longer timcscalcs. In fact, this 
field affects where the base of the convection zone is lo- 
cated. The temporal characteristics of the deeper field is 
still unknown, since the radius data currently available 
have a scant number of data points, poorly distributed 
in time. The amplitude of the variations of the solar 
luminosity produced by the shallower field is of the or- 
der of 0.1%, while the deeper field can produce changes 
of the same order of magnitude. The solar luminosity 
when both magnetic components are strong would be 
about 0.2% higher than when both are weak. 

For instance, the luminosity deficit caused by the 
shallow field at the time of the Maunder minimum is 
about 0.1%. To cause the decrease of the radius ob- 
served between 1715 and 1979, the magnetic field at the 
base of the convection zone had to decrease from about 
1.3 MG to about 600 kG. Since 1979 is at the maxi- 
mum of Cycle 21, the corresponding magnetic field at 
the minimum is smaller than one half of the field at the 
maximum. We thus infer that the magnetic field at the 
minimum is smaller than 300 kG, which is in agreement 
with helioseismology [Antia et al., 2000]. The luminos- 
ity due to the deep field in 1715 is higher than that in 
1979 by about 0.1%, but is almost the same as that in 
1925. The Maunder minimum terminated around 1715, 
when both shallow and deep magnetic field components 
are, according to our model results, strong. We know 
that the shallow magnetic field was weak during the 
Maunder minimum, because of the small sunspot num- 
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Figure 4. Solar variability corresponding to the prob- 
able change in the solar radius between 1715 to 1979. 
Rcz is the location of the base of the convection zone, 
while Rb is the location of maximal magnetic field. MG 
stands for Million Gauss. The mass depth is defined as 
log(l — M/Mq) by the mass coordinate M. The smaller 
the mass depth, the closer to the surface of the model 
sun. 

bers during that period. If we assume that the deeper 
component was also weak, the solar luminosity during 
the Maunder minimum may have been about 0.2% lower 
than that after 1715. Such a large irradiancc change 
may account for the transition between the little ice 
age and the following warm period. 

The solar radius change found to occur between 1925 
and 1976 (or 1979) brings up an additional issue. The 
implication is that the solar luminosity in 1925 was 
about 0.1% higher than that in 1976 (or 1979). Of 
course, the climate consequences of this higher lumi- 
nosity period depend strongly on the duration of this 
phenomenon, which at present is unknown. However, 
there is a possibility that the solar effect alone would 
have caused some cooling, and that the global warming 
caused by the increase of the concentration of green- 
house gases could be more significant than what has 
been estimated ignoring the solar effect. Because of the 
complexity of the climate response to solar irradiance, 
and in view of the serious uncertainties in the current so- 
lar data base, we cannot, at present, draw any sweeping 



conclusions. What we can conclude with some certainty 
is that solar variability could be a significant driver of 
climate change, and that the efforts made to date in ad- 
dressing this issue are inadequate. The physical model 
that we have developed, with all its inadequacies (for 
example, it is unidimcsional, whereas the problem is 
clearly three-dimensional) can constitute a framework 
which could form the basis for interpreting a variety of 
historical solar observations. 
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